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Using mean-field theory, we study the equilibrium properties of boson-fermion mixtures confined 
in a harmonic pancake-shaped trap at zero temperature. When the modulus of the s-wave scatter- 
ing lengths are comparable to the mixture thickness, two-dimensional scattering events introduce a 
logarithmic dependence on density in the coupling constants, greatly modifying the density profiles 
themselves. We show that for the case of a negative boson-fermion three-dimensional s-wave scat- 
tering length, the dimensional crossover stabilizes the mixture against collapse and drives it towards 
spatial demixing. 
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I. INTRODUCTION 



Fermionic atomic gases were brought together with bosonic atoms to quantum degeneracy in several alkali atom 
mixtures, such as 7 Li- 6 Li [l J3] > 23 Na- 6 Li 0], 87 Rb- 40 K 0, [H, and very recently in a mixed gas of ytterbium 
(Yb) isotopes, 174 Yb- 173 Yb Q. The boson-fermion (BF) coupling strongly affects the equilibrium properties of the 
mixture and can drive quantum phase transitions, as collapse [6| in the presence of attractive BF interaction, or spatial 
demixing as recently observed in the context of three-dimensional (3D) atomic fermion - molecular boson mixtures 
[1, [t| , where the strong interspecies repulsion leads to phase separation. 

Such mixtures can be realized from an imbalanced two-component Fermi gas ( 40 K- 40 K or 6 Li- 6 Li mixtures) where 
all minority fermions become bound as bosons and form a Bose-Einstein condensate (BEC). Though imbalanced 
Fermi gases allowed to observe spatial phase separation between bosonic dimers and fermions, the advantage of a two 
atomic species BF mixture is that boson-boson (BB) and BF interactions can be driven independently and that one 
can access attractive BF interactions [13, [H| • 

The structure and the stabilityof trapped BF mixtures were studied in 3D by using the Thomas-Fermi approxima- 
tion for the bosonic component [13 , Il3j and by using a modified Gross-Pitaevskii (GP) equation for the bosons which 
self-consistently includes the mean- field interaction generated by the fermionic cloud [l5[ . Effects of the geometry 
induced by the trap deformation were studied in the Thomas-Fermi regime in a quasi-3D limit, i.e. when collisions 
can be still considered as three-dimensional (l6| . Such a simple model predicts, in a pancake-shaped trap, that the 
stability of the mixture depends only on the scattering length and the transverse width of the cloud. One should 
expect, in a true dimensional crossover, namely including dimensional effects in scattering events, that the mixture 
stability depends critically on the energy, and thus on the number of particles. 

The dimensional crossover from a 3D to a 2D trapped mixture may be studied in the experiments by flattening 
magnetic or dipolar confinements 



17 , or by trapping atoms in specially designed pancake potentials, as rotating traps 
19], rf- induced two-dimensional traps [2Q| or in one-dimensional lattices [2l| where 



, gravito-optical surface traps 
a 3D gas can be split in several independent disks. 

In the limit where scattering events are bi-dimensional, it is well known that a hard-core boson gas shows very 
different features from its 3D counterpart. In 3D, particle interactions can be described by the zero- momentum and 
zero-energy limit of T- matrix, leading to a constant coupling parameter. In 2D, T-matrix vanishes at low momentum 
and energy [22l . [231 ] and the first-order contribution to the coupling is obtained by taking into account the many-body 
shift in the effective collision energy of two-condensate atoms (24|, [25| . This leads to an e nerg y dependent coupling 
parameter that greatly affects the equilibrium and the dynamical properties of the gas [26|, |27| . 

In this paper we study the equilibrium properties of a mixture of condensed bosons and spin-polarized fermions, 
through the dimensional crossover from three to two dimensions, by following the procedure outlined by Roth [3] for 
the 3D mixture. We neglect fermion-fermion interactions and we include BF s-wave interaction self-consistently in 
a suitably modified GP equation for the bosons. For the case of BF repulsive interaction, the increasing anisotropy 
softens the repulsion, and a quasi-3D spatially demixed mixture is mixed in quasi-2D. For the case of BF attractive 
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interactions, the dimensional crossover acts as a Feshbach resonance and induce repulsive interactions, so that a Q3D 
mixture near collapse can be driven towards spatial demixing in Q2D. In the strictly 2D regime the results depend 
on the model one assumes for the bi-dimensional scattering lengths. 

The paper is organized as follows. In Sec. |TT] we introduce the theoretical mean-field model for the description 
of ground-state density profiles of the BF mixture. The models for the coupling through the dimensional crossover 
are outlined in Sec. IIIII The density profiles obtained for a 6 Li- 7 Li and a 40 K- 87 Rb mixtures are shown in Sec. IIVI 
Section fVl offers a summary and some concluding remarks. 

II. MEAN-FIELD MODEL FOR THE DENSITY PROFILES 

We consider a BF mixture in a 2D geometry, with respective particle numbers Nb and Np, confined in harmonic 
trap potentials Vb,f = ,f^b f 7 * 2 - Here tyib,f is boson (fermion) mass and ujb.f is the radial trap frequency as 
seen by boson or fermion species. Within the mean-field approach the total energy functional at T = is written as 

E[^ B ,M = [ d 2 r^-\V^ B ? + V B {r)\^ B \ 2 + \ 9 bbWb 

d 2 r{Tp + V F (r)\yj F \ 2 } + [ d 2 rg BF \M 2 \^F 



(1) 



where i\>b,f is the ground-state wave function of bosons and fermions, respectively. In the above boson species are in 
the condensed state and fermion species is assumed to be spin-polarized and its kinetic energy is written within the 
Thomas- Fermi- Weizsacker approximation as (2|| [2j| 

h 2 ( 2 X w \Vn F \ 2 \ 

T F = nn F + — , (2) 

m F \ is np J 

where np = |V>f| 2 is the fermion density and the Weizsacker constant is Xw — 1/4. Normalization conditions for 
Nb bosons and Np fermions read Nb = J d 2 r\ipB\ 2 and Np = f d 2 r\ipp\ 2 . The interaction couplings between the 
bosons and between bosons and fermions are denoted by gsB and gBF, respectively. One notable difference between 
the form of the energy functional given above and that in 3D, is that the BB and BF interaction strengths are in 
general density dependent in contrast to the situation in 3D. More specifically, in 3D the interaction strengths are 
proportional to the scattering lengths o.bb and asF whereas in 2D as we shall explain below they depend on the 
density or equivalently the chemical potential. The Euler-Lagrange equations for the mixture read [291 ] 

h 2 1 

V 2 + V B + g BB \^B\ 2 + 9bf\^f\ 2 - Ms >i/>b = , (3) 



2ms 



and 

H 2 h 2 1 

•- A W V 2 + V F + 2^ F | 2 + g BF \^b\ 2 - MF # = , (4) 

Imp rap J 

in which we have introduced the chemical potentials (J,b,f for bosons and fermions. The above equations of motion 
are obtained by functional differentiation from E[ips, V'f] neglecting the higher-order terms involving 5g / 8^jb,f which 
is valid in the dilute gas limit tibObb ^ 1 an< ^ n FO- B F *C 1. The dilute gas conditions above further maintain 
that beyond mean- field corrections are not called for. They can become notable when asB, o-bf and/or Nb, Np are 
large for fixed trap frequencies. For the systems under consideration we have chosen the parameters appropriately 
and verified by numerical calculations so that ubo,bb^ n F a BF ^ 1- Therefore, in the examples we shall discuss 
subsequently, the beyond mean-field terms in the energy functional are not important. 

It should also be noted that the existence of BEC in 2D needs to be treated carefully. Initial attempts have 
concluded that no BEC could occur in 2D trapped gases but recent considerations within the Hartree-Fock-Bogoliubov 
approximation, [3(| the density dependent interaction strength [3lJ and numerical simulations [13] have established 
firmly the occurrence of BEC for such systems. Thus, our assumption of a 2D condensate at T = is justified. 



III. 2D INTERACTION MODELS 



In cold atom experiments a 2D geometry is obtained by trapping the atoms in a highly anisotropic trap where the 
axial confinement is very tight, so that the axial potential is on the same order or larger than the chemical potentials 
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of the two components. Within this condition, the axial widths are on the order of the oscillator lengths for the axial 
direction a,j Z = yJh/mjUJjz, tOj Z being the axial trap frequency for bosons (j = B) and for fermions (j = F). For 
simplicity, here and in the following we assume that as z = &Fz = o, z . 

The value of a z with respect to the modulus of the 3D scattering lengths, determines whether the scattering events 
occur in 3D or in 2D, and thus suggests how to calculate the many-body interaction potentials. 

The interaction couplings gsB and gsF are determined microscopically from the effective interaction potentials 
(two-body scattering amplitude, T-matrix, etc.) in the limit of low energy and momenta. In the case of a 3D system, 
the scattering amplitude and qbb and gsF are constants determined by the s-wave scattering lengths cibb and cibf- 
In 2D the scattering theory approaches give rise to a logarithmic dependence [22l . [23| . Starting from a 3D system 
and increasing the anisotropy (by increasing the trap frequencies in the axial direction) the geometry flattens to take 
a pancake shape and eventually a genuine 2D system is obtained. In the following we identify different scattering 
regimes depending on the relation between the axial confinement length and scattering lengths and provide expressions 
for the interaction couplings in these regimes. 



A. Quasi-3D scattering 

In this regime, the axial oscillator length a z of the mixture is assumed to be larger than the modulus of ass and 
clef, the s-wave scattering lengths for BB and BF interactions, respectively. 

The effective BB interaction strength can be obtained by multiplying the 3D value of the coupling g B D B = Anh 2 a/mB 
with a factor |0(O)| 2 = l/-\/27ra z , 4>{z) being the axial wavefunction. This is obtained by assuming that the motion 
in the z-direction is frozen in the ground state of the harmonic potential with trapping frequency u>j Z and integrating 
the 3D GP equation over z (after multiplying with <fc*(z) in the spirit of taking an expectation value. The chemical 
potential /i gets shifted by hwj Z /2). Assuming that the profile for fermions to also be Gaussian in the z-direction, we 
apply the same idea to the BF interaction g B D F — 2nh 2 a/m B p, where tubf is the reduced mass. Thus, we obtain 

2V2nh 2 a B B , y/^h 2 a B F ftL , 

g B B = , and g B F = , (5) 

rriB a z tubf a z 

as the effective interaction couplings in the quasi-3D scattering regime. 



B. Strictly-2D scattering 

This regime cor resp onds to the limit a z <C |c&bb|> \o-bf {■ The coupling parameter we use is from a T-matrix 
calculation [24], [25|, |33| which takes into account the many-body shift in the effective collision energy of two condensate 
atoms and it becomes a self-consistent problem. Since a z < \a BB \, Wbf\ the calculation is purely 2D, and the 



interaction strengths do not depend on the parameters in the z-direction. Al Khawaja et al. 24j argue that when two 
condensate atoms collide at zero momentum they both require an energy \ib to be excited from the condensate and thus 
the many-body coupling is given by evaluating at — 2/j,b the two-body T-matrix (T^) setting gsB = (0|T2{,(— 2/ib)|0). 
On the other hand, Lee et al. [33j calculate T2b at —/is arguing that this result includes the effect of quasiparticle 
energy spectrum of the intermediate states in the collision. Gies et al. (34[ claim that Al Khawaja et al. [2J| argument 
that the excitation of a single condensate atom is associated with an energy of —[Lb includes on the mean-field energy 
of initial and final states and neglects the other many-body effects on the collision which presumably are included in 
the result gBB = (0|T 2 b(— ^s)|0). With this proviso we take 



and similarly, 



9 BB - w 2 77fc2\ ( b ' 

m B in(nBm B a BB /4:h z ) 



9BF = ~ 2lrh2 I ( 7 ) 

rriBF hi((^B + fjLF)m B Fa 2 BF /^h 2 } 



where the scattering lengths asB = Q?bb an< ^ aBF = a BF are m P rm ciple 2D scattering lengths. Our choice for the 



2D scattering lengths will be discussed in Sec. HVl In the case of BF interaction strength we used the reduced mass 
ttibf and made the replacement /i — *• (/is + /^i?)/2. Similar considerations to write down the BF T-matrix were also 
made by Mur-Petit et al. [Hj]. 
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In this regime the interaction parameters must be determined self-consistently. The way the equations are written 
above suggests to solve for wave functions for given values of qbb and gBF, then calculate chemical potentials 



and 



»f = TT I d 2 r[ 7 ^—X w \^ F \ 2 + l:rnFU JF r 2 \ijF\ 2 + — 2^ F \ A + 9bf\^f\ 2 \^b\ 2 \ (9) 

IN F J y 2m p 2 VflF J 

and check whether the expressions for gsB and gsF are satisfied. To follow common practice we start with initial 
chemical potentials, calculate g's and then calculate chemical potentials using the obtained wave functions and require 
self-consistency. Note that in this regime the results do not depend on the value of u>j Z , i.e. on the value of the 
anisotropy parameter A = lobz/^b- 



C. Quasi-2D scattering 

When a z > |<Xbb|) |qsb.f| collisions are bi-dimensional but influenced by the z-direction. In this regime, which is 
in between the previous cases, the 2D scattering length can be expressed in terms of the 3D scattering length [3(|. 
Substituting 



= 2V2J-a z e (10) 
V & 

in the coupling strength expressions for strictly 2D regime with B as 0.915, the coupling strengths now become [33| 

2y / 27r/r ! a BB 
itlb a*z 

1 + -j= Ofs- ln{Bh 2 /2irnBm B a 2 g ) 

and 



9bb - — — i &bb r^ fc 9 Q /o_.. ~ —\ ( n ) 



9BF - — 1 a Ft f 1 / t~> fc 9 /o / . . .. loT' y 11 ) 



IV. RESULTS AND DISCUSSION 



We first consider a lithium mixture with particle numbers Nb = 10 6 and Np — 5 x 10 5 , and radial trapping 
frequencies ujbI^-k = 4000Hz and uif/2tt = 3520Hz. The BB and BF scattering lengths are taken as agg = 5.1 do 
and cibf = 38 clq, respectively, in which ao is the Bohr radius. 

In Fig. [T] we show the density distributions ub{t) and njr(r) of bosonic and fermionic components in the three 
scattering regimes: the quasi-3D, where the coupling is given by Eq. ([5]), the quasi-2D, where the coupling is given 
in Eqs. (jTTJ) and (|12|). and the strictly 2D, where we use the coupling given in Eqs. ([6|) and ((7]) and where we set 
the bi-dimensional scattering lengths equal to af® (Eq. (|10p ) evaluated in the limit of vanishing a z /afP. This choice 
assures the strictly 2D model to be the limiting case of the Q2D, that depicts the crossover behavior. 

When aBF/ciz = 0.1 (top panel) the mixture has 3D character in terms of collisions even though the geometri- 
cal confinement (A = 10 3 ) renders the system 2D kinematically. The calculated chemical potentials Hb/^Bz and 
fip/f^Fz being less than unity also confirms that the system is geometrically 2D. In this regime the density distribu- 
tions for quasi-3D and quasi-2D models look very similar. The boson and fermion components occupy the inner and 
outer parts of the disk giving a segregated phase for the chosen parameters. The 2D model is evidently inapplicable 
in this regime because a,BF/a z < 1. 

In the middle panels of Fig. [T] we show density profiles for the same mixture with asF/a-z = 1 for an anisotropy 
parameter A = 10 5 . This corresponds to a completely frozen motion in the z-direction and to the crossover in the 
scattering properties from 3D to 2D. Figures QJc) andQJd) reveal that the density profiles in the three models are 
very similar, except for the fact that the 2D model predicts a larger spatial extension of the density profiles. 
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aBF/a z = 0.1 





r/a B ± r/a B ± 

FIG. 1: (Color online) Boson and fermion density profiles for 6 Li- 7 Li mixture with Nb = 10 6 and Nf — 5 x 10 5 , radial trapping 
frequencies uib/^tv — 4000 Hz, luf/^h — 3520Hz and scattering lengths ass = 5.1ao, clbf = 38ao where ao is the Bohr 
radius. The length unit is the radial harmonic oscillator length for bosons aai = sfh/mBWB- The density given is in units of 
10~' l aZ 2 L and is normalized to unity. The three regimes agf /a z = 0.1, 1, 10 correspond to values of the asymmetry parameter 
A « 10 3 , 10 s , and 10 7 , respectively. 
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FIG. 2: (Color online) Same as in Fig.Q]for a 40 K- 87 Rb mixture with N B = 10 6 , N F = 5 x 10 5 . The values of \a BF \/a z = 0.3, 10 
correspond to A tti 2 x 10 2 and 2 x 10 5 . 



Finally, in the bottom panel of Fig.[T]we consider aBF/az — 10 with A = 10 7 . a z being smaller than in the previous 
case, the bi-dimensional scattering lengths are smaller and both the 2D and Q2D models predict a mixed phase even in 
the center of the trap, while the Q3D curves still show phase separation. For this anisotropy parameter, the scattering 
events should be truly 2D and our corresponding model should yield the most accurate density profiles. Evidently 
the Q3D model is not yet valid, but we plot it just to compare the predictions of the different models. 

We now turn our attention to 40 K- 87 Rb mixture having an attractive BF scattering length. We consider a system 
with particle numbers = 10 6 and Np — 5 x 10 5 , and radial trapping frequencies lob/2-k = 257Hz and uif/2tt = 
378 Hz. The BB and BF scattering lengths are taken as clbb = HO do and clef — — 284 do, respectively }37j . 
For attractive interactions, the effective 2D BF scattering length is positive [see Eq. (fTTJ|) ]. namely the dimensional 
crossover induces effective repulsive interactions (35j | , as alr eady predicted in a condensate with attractive boson-boson 
interaction 36]. Thus, the strictly 2D couplings (Refs. [24L l33l [34 [3oT|) refer to hard-core collisions (22I ]. 

Figure [5] illustrates the density profiles nsir) and np(r) in quasi-3D and 2D scattering regimes, characterized by 
o-bf/o-z = —0.3 (A k 2 x 10 2 ) and aBF/a z = —10 (A m 2 x 10 5 ), respectively. In the case asF/a z = —0.3, we 
observe that the density profiles are similar for quasi-3D and quasi-2D models and show a bump in the center of 
the fermionic density due to the attractions with the bosons. For asF/o-z = —10, the Q2D model approaches the 
2D one, the only difference being that the first model predicts complete spatial separation between the bosonic and 
the fermionic components, while the second predicts a residual mixed phase at the center of the trap. The crossover 
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between the two regimes is shown in Fig. [3] For A < 10 2 , the fermionic density is enhanced at the center of the trap 
because of the presence of the bosons. In this regime the BF coupling term is negative, as shown in Fig.HJ At A = 10 s 
the fermions are pushed out of the center of the trap because of the large repulsive BF interaction (see Fig.[4j. By 
increasing further and further the anisotropy, the BF coupling is still positive but decreases and the two components 
are partially mixed. For 10 2 < A < 10 5 no stable solutions are found. Thus, as shown in Fig.[U the dimensional 
crossover plays the role of a Feshbach resonance. Squeezing the trap one may naively expect the gas just collapsing, 
but the crossover in the scattering geometry changes the nature of the instability from collapse to demixing, and a 
further squeezing of the trap stabilizes the mixture. All curves shown in this sections correspond to densities that 
fulfill the diluteness conditions ubObb ^ 1 an d nFa 2 BF <C 1, even at close to the resonance shown in Fig.H 



V. SUMMARY 



In summary we have studied the equilibrium properties of a boson-fermion mixture confined in a pancake-shaped 
trap, in the dimensional crossover from 3D to 2D. The boson-boson and the boson-fermion couplings used are those 
derived from the two-body T-matrix evaluated (i) at zero energy in 3D, (ii) taking into account the discreteness of the 
spectrum in the axial direction, in the crossover, (iii) taking into account the many-body energy shift in the strictly 
2D limit. The density profiles and the couplings have been evaluated self-consistently using suitable modified coupled 
Gross-Pitaevskii equations for the bosonic and the fermionic wave functions. 

For the case of a positive 3D boson-fermion scattering length, the dimensional crossover softens the repulsion, 
so that the components of a demixed boson-fermion mixture in 3D can mix in the 2D limit. For the case of a 
negative 3D boson-fermion scattering length, the dimensional crossover is more dramatic and plays the role of a 
Feshbach resonance. Our study shows that the squeezing of the pancake-shaped trap may drive a strong-attractive 
unstable mixture towards a stable mixed mixture passing through a demixed phase. This numerical study may be 
reproduced in the actual experiments with BF mixtures. The goal being to reach a regime where the modulus of 
the scattering lengths is comparable or greater than the mixture axial size, one may exploit Feschbach resonances 
to increase the magnitude of the 3D scattering lengths, or one may engineer very flat traps as already done in the 
context of experiments with a single BEC component. 
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FIG. 4: (Color online) Effective BF interaction strength (in dimensionless units) for the 40 K- 87 Rb mixture within the Q2D 
model as a function of the anisotropy parameter A. Dots refer to the numerical calculation performed in the Thomas-Fermi 
(TF) approximation, namely neglecting the Laplacian terms in Eqs. Q and (UJ), while triangles refer to the full solution of the 
same equations (GPE). 
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